The /-electron challenge: localized and itinerant states in lanthanide oxides united by GW@IA}A+U 
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Many-body perturbation theory in the GW approach is applied to lanthanide oxides, using the local-density 
approximation plus a Hubbard U correction (LDA-l-i7) as the starting point. Good agreement between the Go Wo 
density of states and experimental spectra is observed for Ce02 and Ce203. Unlike the VDA+U method Go Wo 
exhibits only a weak dependence on 17 in a physically meaningful range of U values. For the whole lanthanide 
sesquioxide (Ln203) series Go Wo@LDA-l-(7reproduces the main features found for the optical experimental 
band gaps. The relative positions of the occupied and unoccupied /-states predicted by Go Wo confirm the 
experimental conjecture derived from phenomenological arguments. 

PACS numbers: 71.10.-w,71.15.-m,71.20.-b,71.27.+a 



The accurate first-principles description of the electronic 
structure of /-electron systems, i.e. materials containing lan- 
thanide or actinide elements, is currently regarded as one of 
the great challenges in condensed matter physics, /-electron 
systems are characterized by the simultaneous presence of 
itinerant (delocalized) and highly localized /-states and in- 
teractions between them. Most computational methods are 
suited only to one type. Density-functional theory (DFT) - 
currently the standard approach for electronic structure calcu- 
lations of extended systems - proves to be inadequate for /- 
electron systems in the most commonly applied local-density 
or generalized gradient approximation (LDA or GGA, respec- 
tively). One of the major deficiencies of LDA and GGA is the 
delocalization (or self-interaction) error which is particu- 
larly severe for systems with partially occupied d or /-states 
and can even lead to qualitatively incorrect metallic ground 
states for many insulating systems. Hybrid functionals yfl, on 
the other hand, partly correct the self-interaction error by in- 
corporating a certain portion of exact exchange, which signif- 
icantly improves the descriptions of d- or /-electron systems 
Jsl Q]- The dependence on adjustable parameters, however, 
remains a concern. Conversely, correlation effects that govern 
the physics of localized /-electrons can in principle be treated 
systematically by dynamical mean field theory (DMFT) |5]. 
In practice these many-body coiTections are only applied lo- 
cally to an atomic site (e.g. the Anderson impurity model) 
and the impurity solvers require input parameters (such as the 
Hubbard U) for the interaction strength. Moreover, most ex- 
isting DMFT schemes are coupled (non self-consistently) to 
local or semilocal DFT calculations and the description of the 
itinerant electrons therefore remains on the level of LDA and 
GGA. 

As a first step towards a systematic ab initio understand- 
ing of /-electron systems, we apply many-body perturbation 
theory (MBPT) in the GW approach to a selected set of lan- 
thanide oxides (Ce02 and Ln203 (Ln=lanthanide series)) |Q] 
in this Letter. These compounds have important technolog- 
ical applications 10. S S U. 12], in particular in catal- 
ysis, where Ce02 -based compounds have attracted consider- 
able interest from both experiment and theory I sl [Tsi \\A [Tsi 



ilESIilHIIS]- Unlike in most previous studies, the GW 
calculations in this work are based on LDA ground state cal- 
culations including a Hubbard U correction (henceforth de- 
noted GoWo@LDAh-C/). Our GoWo@UiK+U calculations 
provide a qualitative understanding of the general trend ob- 
served for the band gaps of the Ln203 series and reproduce 
the characteristic features of the series, in particular the four 
dips observed in the experimental curve. 

The GW approach corresponds to the first order term of 
a systematic expansion in MBPT [21J and has become the 
method of choice for the description of quasiparticle band 
structures in weakly correlated solids [22]. Through the 
screened Coulomb interaction W it captures the screening 
among itinerant electrons while at the same time treating ex- 
change at the exact exchange level (given by the Hartree-Fock 
expression). The latter should account for a large part of the 
many-body interactions among localized d or /-electrons, as 
demonstrated recently for CU2O and VO2 ll23ll . In this Let- 
ter we challenge the conventional view that regards many 
/-electron systems as strongly correlated electron systems 
for which band theory is inadequate [24]. In the hierarchy 
of many-body perturbation theory, strong correlation denotes 
correlation effects that go beyond exact exchange and the 
weak coiTelation regime of GW . The good agreement be- 
tween our GW calculations and available experimental data 
demonstrates that the GW method can treat both itinerant spd 
bands and localized /-bands accurately for the materials we 
have considered. Since the screened exchange picture of the 
GW approach captures the essential physics we challenge the 
classification of these materials as strongly correlated. 

The GW method is typically applied in a perturbative man- 
ner (henceforth denoted GqWq) in which the quasiparticle 
(QP) energies are calculated as a first-order correction 
to the eigenenergies e„k and eigenvectors V'nk of a reference 
single particle Hamiltonian Hq 



V^^ I^Ank) (1) 



Here YT^'^ is the GqWq self-energy calculated from the one- 
particle Green's function Go and screened Coulomb interac- 
tion Wo, both evaluated using e„k and i/'nk, and V^'^ is the 
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FIG. 1: (Color online) The DOS of Ce02 (left) and CeaOa (right) from LDA, LDA+U and GoH^o@LDA+;7 with U=5A eV are compared to 
experimental data(XPS+BIS and XPS+XAS) t25J. The LDA and LDA+17 curves have been offset vertically for clarity. 



exchange-correlation potential included in Hq- For most of 
the sesquioxides considered in this work LDA and GGA in- 
correctly predict a metallic ground state. In these cases first 
order perturbation theory based on LDA or GGA is not ap- 
plicable and alternative reference Hq have to be employed 
Sll^llilli]. In this work we use the LDA+U method |30] 
as the starting point for GqWq. By adding a site- and orbital- 
dependent correction 6V^ to the standard LDA single-particle 
Hamiltonian, LDA+U significantly improves the description 
of highly localized states, and therefore overcomes the major 
failure of LDA for these systems. To describe highly local- 
ized states accurately, we have implemented an all-electron 
GqWo approach 13 If based on the full-potential linearized 
augmented plane wave method ll32ll . 

The LDA+U method is conceptually similar to GW. It is, 
however, not a substitute for GW even for localized states: 
1) The link between the LDA+U and the GW approximation 
relies on the assumption that the hybridization between lo- 
calized and itinerant states can be neglected, which in many 
cases is not valid; 2) the 6V^ correction in LDA+U has di- 
rect effects only on the corresponding localized states; the 
description of itinerant states remains at the LDA level; and 
3) screening in LDA+U is static, while in reality screening 
is dynamic and, has a stronger energy-dependence for local- 
ized electrons than for itinerant ones. The LDA+U approach 
by itself is therefore not expected to provide a satisfactory de- 
scription to the electronic structure of /-electron systems. 

An advantage of the GoWo@LDA+U approach lies in the 
fact that the Hubbard U corrections enter self-consistently in 
the ground state calculation. This becomes important when 
localized states hybridize with band states. A less appealing 
aspect of the LDA+U approach concerns the parameter U, 
which, in many cases, is determined by fitting to experimental 
data. The onsite Coulomb interaction U, however, has a well- 
defined physical meaning, and can be calculated from first- 
principles (see, e.g. Refs. Ii33il ). In addition we demonstrate 
below that GqWo based on LDA+U is much less sensitive 
to U than the LDA+U itself. The GqWq calculations also 
remove the problem of the double counting corrections that 



are not well defined in the LDAH-J7 approach, by subtracting 
them out as part of Vxc in Eq. [T] Despite these advantages the 
GqWo @LDAh-J7 approach will not be suitable in cases where 
strong correlation effects become important (e.g. the Kondo 
resonance), for which many-body interactions that go beyond 
the GW approach have to be included. This is currently the 
domain of DMFT, as alluded to in the introduction. 

Figure [T] shows the density of states (DOS) of Ce02 and 
CesOa calculated from LDA, LDA+U, and GoWo@LDA+U 
(with U = 5.4 eV ll34ll ) together with the experimental spec- 
tra. The GqWq density of states for Ce02 agrees well with 
the experimental data from direct (XPS) and inverse (BIS) 
photoemission spectroscopy or X-ray absorption spectroscopy 
(XAS). In Ce02, the empty /-states introduce a sharp peak 
in the fundamental band gap formed between the 0-2p va- 
lence and Ce-5rf conduction band. The most important quan- 
tities here are therefore the p-f and p-d gaps. For the latter 
the GqWq value of 6.1 eV is in good agreement with the ex- 
perimental one of 6.0 eV. The p-f gap, however, cannot be 
unambiguously determined from XPS-BIS or other available 
measurements (See, e.g. Ref. |18] and references therein). 
As expected, LDA underestimates both gaps, but the p-d gap 
is only slightly smaller than in experiment (5.5 vs 6.0 eV), 
whereas LDA+U decreases it to 5.1 eV. We also note that our 
GqWo @LDA+U DOS for Ce02 is very close to that obtained 
from the recently proposed self-consistent GW method I35I1 . 

More intriguing features are observed for Ce203. As ex- 
pected, the on-site Hubbard correction in the LDA+U splits 
the single /-peak in LDA to occupied and unoccupied /- 
bands (denoted as f°'^'^ and /"" , respectively), the former 
falling within the p-d gap and the latter overlapping with the 
conduction bands. The p-d gap is nearly the same as in LDA, 
but the p-f°'^^ splitting is greatly reduced. Applying the 
GqWq corrections to the LDA+U ground state, we observe 
two remarkable features: 1) the 0-2p band is pushed to lower 
energy with respect to the f°^'^ band, and 2) the /"" -band 
shifts up in energy away from the conduction band edge in- 
creasing the splitting between the J°'^'^ and /"" bands at the 
same time. 
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FIG. 2: (Color online) Band gaps of Ce02 (p-/ and p-d, up- 
per row) and Ce203 (/-d and p-d, lower row) from LDA-l-f/ and 
GoW^o@LDA-l-?7 as a function of (7. The dashed rectangles in the 
lower panel indicate the range of physically meaningful values of \J 
(see text). 



The GoVFo and the experimental spectrum in Fig. [T]are 
aligned at the upper valence band edge and not the Fermi level 
since this is not well defined at the K at which the calcu- 
lations are performed. We again find the peak to be in 
good agreement with the BIS data. The position of the 
peak, however, differs by approximately 1 eV. With respect to 
the band gap the comparison with experiment is aggravated 
by the limited experimental resolution. Taking the difference 
between the conduction band edge and the upper edge of the 
XPS-XAS peak gives a band gap that is consistent with 
the optical band gap of ~2.4 eV [St |36|l and the GqWq gap 
of 2.0 eV. If one instead references against the peak center of 
the states, the experimental band gap of Ce203 would 
be larger than 3 eV. Further experimental evidence is clearly 
needed to determine the actual value. 

Figure |2] illustrates the influence of U on the LDAh-L/ and 
GoWo@LDAh-[/ calculations for the examples of Ce02 and 
Ce203. Since the /-states are essentially empty in Ce02, 
the effect of U is relatively weak: For C/=0 to 8 eV, the p- 
/ and -p-d gaps from LDAh-[/ change only by approximately 
1.0 and 0.5 eV, respectively, and become nearly constant in 
GoW^o- The situation is more complex in Ce203. In LDAh-C/, 
the /-d gap depends sensitively on U and varies by nearly 3 
eV, but the -p-d gap remains almost unaffected. In contrast, 
both the j-d and the p-d gap exhibit a slight U -dependence 
in GqWq, changing by 1.2 and 0.6 eV over the full U range 
explored here. Most importantly, however, the U -dependence 
reduces to only approximately 0.3 eV in the range of "phys- 
ically meaning ful" values of U (-5-7 eV ilMlli), which is 
akeady in the range of experimental error bars. We note in 
passing that the relatively weak dependence on U has also re- 
cently been reported by Kioupakis etal., who applied a similar 
approach to solid hydrogen [37]. In addition. Figure |2] shows 
that within LDAh-[/ one could obtain an apparently more ac- 
curate p-/ (in Ce02) and j-d (in Ce203) gap by using a sig- 
nificantly larger L/, which, however, would not improve the 
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FIG. 3: (Color online) Band gaps of the Ln203 series from hDA+U 
and Go Wo ((7=5.4 eV) are compared to SIC-LDA results 1 10] and 
experimental optical gaps [8]. The schematic in the upper part of the 
figure illustrates the position of the f°'^'^ and states extracted 
from the Go Wo QPhDA+U calculations in relation to the valence and 
conduction band edge (VB and CB). 



description of itinerant states including, in particular, the p-d 
gap- 
It has long been recognized that, although many proper- 
ties of rare-earth compounds exhibit a monotonous behavior 
across the lanthanide series, some show a striking variation. 
For example the optical band gaps of rare earth sesquioxides 
[8] shown in Fig. [3] exhibit clear dips for Ce, Eu, Tb, and Yb, 
which appear to be unaffected by structural variation across 
the series. In Fig.[3]the optical gaps of the Ln203 series li36ll 
are compared to LDAh-C/, Gq Wq @LDAh-C/ and previous self- 
interaction corrected LDA (SIC-LDA) results |flot]. Since U 
is the effective e-e interaction among /-electrons screened by 
spd states, which are very similar in all Ln203 compounds, 
we expect U to be only weakly dependent on the number of 
localized /-electrons, and therefore use a constant [7=5.4 eV 
for the whole series. Compounds denoted by filled circles in 
Fig-EJcrystallize preferentially in the hexagonal structure, for 
which all calculations have been performed. Starting from 
Sm203 (denoted by open circles), the most stable phase at 
room temperature is cubic bixbyite, but the middle members 
can also exist in the monoclinic phase [9]. 

As can be seen from Fig. [3] all the essential features in the 
experimental curve are reproduced by the GqWq calculations 
including the four dips and the behavior in between. Even 
the quantitative agreement is good for most compounds. In 
addition our first principles calculations provide easy access 
to the character of each peak in the DOS and thus the char- 
acter of the band gap, which is schematically shown in the 
upper part of Fig. [3] In La203 (empty /-shell) the band gap is 
formed between the 0-2p valence and the La-5d conduction 
band. As the occupation of the /-states increases, both f°'^'^ 
and /"" continuously move downward in energy and the band 
gap evolves from p-d via f-d to p-f. This process repeats 
itself in the second part of the series (starting from Gd203) 
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where the spin-up /-states have become fully occupied (and 
lie deep below the 0-2p states) and the spin-down /-states 
move downward in energy with increasing occupation. 

The character of the band gap across the series agrees well 
with the experimental conjecture derived from phenomeno- 
logical arguments | Tj, |8l (38t] . This is not the case in SIC-LDA 
and LDA+U. Not only does LDA+U underestimate the band 
gaps of most Ln203 compounds, it also only shows a weak 
minimum at Tb203 and fails to reproduce the plateau between 
H02O3, Er203, and Tm203. 

The remaining quantitative differences between the Go Wo 
and experimental curve (especially for the later members of 
the series) could be due to several factors: 1) the experimen- 
tal error bar, sample quality and instrumental resolution are 
not known 181; 2) to exclude the influence of excitonic ef- 
fects the Go Wo results should be compared to photoemission 
(both direct and inverse) spectra, which are unfortunately not 
available for most members of the Ln203 series; 3) spin-orbit 
and multiplet effects are not taken into account in our calcu- 
lations; and 4) our Go Wo calculations are performed for the 
hexagonal structure, but the later members of the Ln203 series 
(after Ln=Sm) crystallize in the monoclinic or cubic bixbyite 
structure [|6fl. Further investigations of these issues will be ad- 
dressed in the future. 

We close this Letter with some remarks about strong corre- 
lation. A system is often commonly (but by no means satis- 
factorily) classified as strongly correlated if DMFT provides 
an improved description over LDA and it is therefore interest- 
ing to compare GoWa@LDA+U and LDAh-DMFT [5]. Both 
methods add dynamic (i.e. energy dependent) effects to the 
LDA+U Hamiltonian. LDAh-DMFT is based on an under- 
lying picture of localized states, in which higher order cor- 
relation effects among the localized states can be easily in- 
corporated. However, delocalized states (e.g. spd states in 
/-electron systems) are still treated at the LDA level. In con- 
trast. Go Wo based on LDA+U improves upon LDA+U in the 
band picture and both localized and itinerant states are treated 
at the same level (i.e. Go Wo). For Ce203 this can be di- 
rectly quantified by comparing the GaWa@LDA+U density 
of states with that from recent LDAh-DMFT calculations by 
Pourovskii et. al. li20ll . The /-states are described in a simi- 
lar manner, but the p-d gap in LDA+DMFT amounts to only 
~3.5 eV, which is significantly smaller than the experimental 
(~ 5.5 eV) and the GoWo@LDA+J7 (- 5.1 eV) value. For 
many open d- or /-shell compounds it can therefore be ex- 
pected that the GW approach can provide a consistent and 
accurate treatment for both localized and itinerant states and 
it will be illuminating to see when many-body effects beyond 
the Go Wo approximation become significant. 
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